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Abstract 

The Peccei-Quinn (PQ) solution of the strong-CP problem requires the existence of axions, 
which are a viable candidate for Dark Matter. Here we show that, if the Nambu-Goldstone 
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{Sj . potential of the PQ model is replaced by a potential V(|3>|) admitting a tracker solution, the scalar 

>. 

field |<3?| can account for Dark Energy, while the phase of <J> yields axion Dark Matter. Such Dark 
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Matter and Dark Energy turn out to be weakly coupled. If V is a SUGRA potential, the model 



the quark-hadron transition, |$| naturally passes through values suitable to solve the strong-CP 



O : essentially depends on a single parameter, the energy scale A. Once we set A * lO^GeV, at 

o 

Q-f problem, later growing to values providing fair amounts of Dark Matter and Dark Energy. In this 

model, the linear growth factor, from recombination to now, is quite close to ACDM. The selected 

■*— > 

A value can be an indication of the scale where the soft breaking of SUSY occurred. 
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PACS numbers: 98.80.-k, 98.80.Cq, 95.35. +d 
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INTRODUCTION 

The solutions of the strong CP problem proposed by Peccei & Quinn in 1977 (|l|, PQ 
hereafter) leads to one of the accepted models of Dark Matter (DM). PQ consider the 
lagrangian term 

Ce = -L6G-G (1) 

(a s : strong coupling constant, G and G: gluon field tensor and its dual), yielding CP viola- 
tions in strong interactions, and show that its effects are suppressed by making 9 a dynamical 
variable, approaching zero in our cosmic era, its residual oscillations appearing as DM J2J |3j . 
The 9 dynamics is set by assuming that a complex field $ = <pe e j\[2 exists, whose 
evolution is ruled by a Nambu-Goldstone potential 

V(|$|) = A[|$| 2 -F^] 2 , (2) 

which is clearly U(l) invariant. At T < Fpq (the PQ energy scale, which shall be ~ 
10 12 GeV), falls into the potential minimum, so that the U(l) symmetry breaks, as 9 
acquires different values in different horizons. When the chiral symmetry is also broken, close 
to the quark-hadron transition, a further term must be added to the effective lagrangian, 
arising because of instanton effects. This term reads 

V x = [J2{0(T)\qq\0(T))m q ] (1 - cos9) . (3) 

q 

At T ~ 0, the square bracket approaches rn%f% (m n , f^: 7r-meson mass, decay constant). 

The choice of a NG potential is the simplest possible. Here we explore the possibility of 

replacing it by a potential with a tracker solution J4J, |^| . Instead of taking a value ~ Fpq 

soon, evolves over cosmological times. As in the PQ case, the potential shall involve 

a complex field $ and be U(l) invariant. While <fi rapidly settles on the tracker solution 

(apart of residual fluctuations) in almost any horizon, the symmetry is broken soon by 

the values taken by 9, which suffers no dynamical constraints and is therefore random, 

in different horizons. Later on, when a mass term arises because of the chiral symmetry 

break, dynamics becomes relevant also for the 9 degree of freedom, as in the PQ case. Here 

this happens while still evolves over cosmological times. Finally, in the present epoch, 

accounts for Dark Energy (DE). Hence, besides of yielding DM through its phase 9, the $ 

field, introduced to solve the strong CP problem, accounts for DE through its modulus 0. 



Within this model, DM and DE will be weakly coupled. If we take a generalization of 
the SUGRA potential [5( as tracking potential, with an energy scale A ~ 10 10 GeV, we find 
reasonable values for today's DM and DE densities, while 9 is driven to values even smaller 
than in the PQ case, so that CP is apparently conserved in strong interactions. In turn, A 
may be an indication of the scale where the soft breaking of super-symmetries occurred. 

LAGRANGIAN THEORY 

The lagrangian C = \/—g{g^ v d l j^d ll ^ — V{\$\)} can be rewritten in terms of (j) and 9, 
adding also the term breaking the U(l) symmetry, as follows: 

£ = V=g{\gAd^<P + <f>%ed v e] - v{<f>) - m 2 (r, 0)0 2 (i - cos 9)} . (4) 

Here g^ is the metric tensor. We shall assume that ds 2 = g^dx^dx" = a 2 (dr 2 — rjijdxidxj), 
so that a is the scale factor, r is the conformal time; greek (latin) indeces run from to 
3 (1 to 3); dots indicate differentiation in respect to r. Around the energy scale Aq C d 
(quark-hadron transition), we shall take p] 

m(T,0)^ 0.1 m o (0)(^) 3 ' 8 (5) 

with m ((f>) = m-xfir/rf)- At T <~ 0.3-0.2 Aqcd, "^{T,^) shall already approach its low-T 
behavior m (<f>). The equations of motion then read 

9 + 2 (- + t) + m 2 a 2 sm9 = , (6) 

+ 2-0 + aV(0) = 0# 2 , (7) 

and will be mostly used with sin# ~ 9. Then, energy densities pg^ = pe^Hn + Pe,<t>;pot an d 
pressures pg^ = pe^kin ~ Pe^pot, under the condition 9 <C 1, are obtainable from 

:, _m 2 (T,0) 2 
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\kin — TT^u ) Pe,pot ~ : V 



2a 2 

(j) 2 

b,kin = Y~2 ' ^>P ot = V\<i>) ■ (8) 



THE CASE OF SUGRA POTENTIAL 

When 9 undergoes many (nearly) harmonic oscillations within a Hubble time, {pe,kin) — 
{Pe, P ot) and (pg) vanishes. Under such condition, using eqs. ©,©,©, it is easy to see that 

(9) 
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When m is given by Eq. (|5|). m/m = —<$/<$ — 3.8 T/T. At T ~ 0, instead, rh/m ~ — <; 
Here below, the indices 9, <p will be replaced by DM, DE. Eqs. (JHJ) clearly show an exchange 
of energy between DM and DE. Let us notice that the former Eq. © can be formally 
integrated, yielding p DM oc m/a 3 . In particular, this law holds at T <C Aq C d, an d then 

Pdmo?4> — const., (10) 

so that the usual behavior Pdm oc a ~ 3 is modified by the energy outflow from DM to DE. 
Let us now assume that the potential reads 

y(0) = ^exp(4vr0 2 /mj) (11) 

and does not depend on 9; in the radiation dominated era, it admits the tracker solution 

r +2 = g a A a+ W , (12) 

with g a = a(a + 2) 2 /4(a + 6). This solution holds until we approach the quark-hadron 
transition. Then, in Eq. (|7jl. the term <j)9 2 due to the DE-DM coupling, exceeds a 2 V 
and we enter a different tracking regime. This is shown in detail in Fig. 1, obtained for 
matter (baryon) density parameters Q m = 0.3 (fl& = 0.03) and h = 0.7 (Hubble constant 
in units of 100 km/s/Mpc). In particular, Fig. 1(a) shows the transition between these 
tracking regimes. Fig. 1(b) then shows the low-z behavior (1 + z = 1/a), since DE density 
exceeds radiation and then gradually overcomes baryons (at z ~ 10) and DM (at z ~ 3). 
Fig. 1(c) is a landscape behavior of all components, down to a = 1. Notice, in particular, 
the a dependence of pdm, occurring according to Eq. (JTUJ. In Fig. 2 we show the related 
behaviors of the density parameters fli (i = r, b, 9, <p, i.e. radiation, baryons, DM, DE). 

In general, once the density parameter Qde (at z = 0) is assigned, a model with dynamical 
(coupled or uncoupled) DE is not yet univocally determined. For instance, the potential (fTTj) 
depends on the parameters a and A and one of them can still be arbitrarily fixed. Other 
potentials show similar features. 
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FIG. 1: Densities of the different components vs. the scale factor a. Fig. 1(a) magnifies the onset 
of the oscillation regime. Fig. 1(b) shows the low— z behavior. Fig. 1(c) is a landscape picture of 
the whole evolution. All abscissas are log(a/oo). 

In the present case such arbitrariness no longer exists. Let us follow the behavior of pdm, 
backwards in time, until the approximation < 1 no longer applies. This moment must 
approximately coincide with the time when 9 enters the oscillation regime. This occurs when 



2(d/o + (j>/4>) ~ m(T, 0) a 



(13) 



(see Eq. |BJ). At that time, according to Eq. (JTUJ), which is marginally valid up to there, and 
taking 6 = 1, 



Pdm — Po,dm 



0(a) a 3 



m 2 [r(a),0(a)]0 2 (a) 



(14) 
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FIG. 2: Density parameters f2 r & # ^ (radiation, baryons, DM, DE) vs. the scale factor a. 

The system made by eqs. ()13|) .([14 |) . owing to Eq. (|12p. yields the scale factor a h when fluc- 
tuations start and the value of A in the potential ()11|). as soon as p ,DM (the present density 
of DM) is assigned. 

The plots shown in the previous section, drawn for CIdm — 0.27, are obtained for A ~ 
1.5- 10 10 GeV, as is required by eqs. (fH?j) . (fb"fj) . In this case a^ ~ 10 -13 . When VL DM goes from 
0.2 to 0.4, log 10 (A/GeV) (almost) linearly runs from 10.05 to 10.39 and ah steadily lays at 
the eve of the quark-hadron transition. A model with DE and DM given by a single complex 
field, based on SUGRA potential, therefore bears a precise prediction on the scale A, for the 
observational Qdm range. In turn, we can say that, if the soft breaking of super-symmetries 
occurred at a scale slightly above 10 10 GeV, Qde ~ 0.3 is a natural consequence. 

EVOLUTION OF INHOMOGENEITIES 



Besides of predicting fair ratios between the world components, a viable model should 
also allow the formation of structures in the world. This matter will be treated in detail 
in a forthcoming paper. Let us however outline that the model treated here belongs to 
the class of coupled DE models treated by Amendola J2| , with a time-dependent coupling. 
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FIG. 3: Time evolution of DM and baryon fluctuations. The top figure shows DM and baryon 
fluctuation evolution in this model. The two bottom figures compare DM and baryon fluctuation 
evolutions in this model (solid curve), ACDM (dot-dashed), coupled DE (dashed). 

In fact, for small $'s, the r.h.s. of eqs. (J7j). after averaging over cosmological times, reads 
C((p)(po)a 2 with C(<f>) = 1/0. Similarly, in eqs. 0, which are already averaged, the r.h.s. 
are ±C(0) cf> pe (C is the DE-DM coupling introduced in |2|). Let us also outline that Fig. 2 
shows a <p-MDE phase, typical of this class of models, after matter-radiation equivalence, 
as the kinetic energy of DE is non-neglegible during the matter-dominated era. 

By solving the fluctuation equations in J2J, with the above C(0), we find the behavior 
shown in Fig. 3 (top). Figs. 3 (bottom) compare fluctuation evolutions in this model (solid 
curves), with those in an analogous ACDM model (dot-dashed curves) and in a coupled DE 
model with constant coupling C = 0.25 \f%ixG ~ (C(0)) (dashed curves). These plots show 
that the linear growth factor, from recombination to now, is significantly smaller than in 
coupled DE models with constant coupling and, more significantly, is quite close to ACDM. 
The essential differences from ACDM are that: (i) objects should form earlier; (ii) baryon 
fluctuations keep below DM fluctuations until very recently. 



DISCUSSION 

The first evidences of DM date 70 years ago, but its non-baryonic nature became com- 
pulsory in the Seventies, when BBNS and CMBR anisotropics were studied. DE is younger, 
but is now required both by SNIa data 8], as well as by CMBR and deep galaxy data 
9j, |l0| • Axions have been candidate DM since the late Seventies, although various studies, 



as well as the occurrence of the SN 1987a, have finally constrained the PQ scale around 
values Fpq ~ 10 12 GeV. Contributions to DM from to pol ogical singularities (cosmic string 
and walls) have also narrowed the constraints to Fpq [11]. Here they were disregarded and 
could cause shifts in quantitative predictions. We shall deepen this point in further work. 

The fact that scalar fields can yield both DM or DE, just changing an exponent in the 
potential, stimulated the work of various authors. A potential like (JTTj) was considered in 
spintessence models |l2[|. According to the choice of parameters, $ yields either DM or DE. 

On the contrary, in this note we deal with the possibility that $ accounts for both DE 
and DM, and that the strong-CP problem is simultaneously solved. As in the PQ model, 
the angle 9 in Eq. (|1J) is turned into a dynamical variable, i.e. into the phase of a complex 
scalar field $, and is gradually driven to approach zero, by our cosmic epoch. Residual 
9 oscillations, yielding axions, account for DM. The critical time for the onset of coherent 
axion oscillations occurs at the eve of the quark-hadron transition, because of the rapid 
increase of m(T, </>). Here replaces the constant Fpq scale. This stage is illustrated by 
Figs. 4(a), 4(b) where the behaviors of p e , po t, Pe,kin, and v<f>,0 = V<t>fil ' P<t>,6, are plotted. 

The novel features of this model arise because the expectation value of </> is not a constant 
F P Q, but evolves over cosmo logical times; in a sense, we let F P q evolve so to yield DE. Such 
evolution modifies the friction term in Eq. (jUJ). The damping of 9 oscillations is therefore 
greater and 9 oscillations are smaller today. Further, accordingly to Eq. (J3J) the axion 
mass (or the oscillation frequency), which varies fairly rapidly during the formation of qq 
condensate, continues to evolve, over cosmological scales, due to the evolution of the <fi field. 
The (low-)z dependence of m o (0) is shown in Fig. 5. We draw the reader's attention on the 
rebounce at z ~ 10, whose implications on halo formation could be critical |l3| |. 

Constraints on PQ axions came from z = observations, which must be fulfilled by the 
same Fpq scale, fulfilling also cosmological requirement. Here, <fi attains values ~ m p today, 
so that most these constraints should be naturally satisfied. This matter, as well as the 
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FIG. 4: The onset of coherent axion oscillations, at the eve of the quark-hadron transition, due 
to the increase of m(T, rf>), causes the behaviors of pg tPO t, Pe,kin ( a ) an d w^,e = Pcp,e/P4>,e (b) shown 
here. 




log(a/a ) 

FIG. 5: (p variations cause a dependence of the effective axion mass on scale factor a, which is 
shown here. 



question of a direct axion detection, will be deepened in further work. 

Let us outline that the choice of a SUGRA potential is arbitrary and could be replaced by 
other potentials, perhaps better approaching data at z — 0. However, using this potential 
allows to appreciate the conceptual echonomy in this approach. In the PQ approach, the 



Fpq scale is assumed. Here, once A is set in a physically significant range, around the 
quark-hadron transition, <fi naturally passes through values enabling to solve the strong-CP 
problem and later naturally grows to values providing fair amounts of DM and DE. 
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